Decay of the Loschmidt echo in a time-dependent environment 
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We study the decay rate of the Loschmidt echo or fidelity in a chaotic system under a time- 
dependent perturbation V(q, t) with typical strength K/tv- The perturbation represents the action of 
an uncontrolled environment interacting with the system, and is characterized by a correlation length 
£o and a correlation time to. For small perturbation strengths or rapid fluctuating perturbations, the 
Loschmidt echo decays exponentially with a rate predicted by the Fermi Golden Rule, 1/f = r c /ry, 
where typically r c ~ mirifro, £a/v] with v the particle velocity. Whenever the rate 1/f is larger than 
the Lyapunov exponent of the system, a perturbation independent Lyapunov decay regime arises. 
We also find that by speeding up the fluctuations (while keeping the perturbation strength fixed) 
the fidelity decay becomes slower, and hence, one can protect the system against decoherence. 

PACS numbers: 03.65.Yz, 05.45.Mt 
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I. INTRODUCTION 

The time-evolution of a quantum system is quite ro- 
bust to changes of the system initial conditions, irrespec- 
tive of the nature of the underlying dynamics |l| . This is 
in deep contrast to classical evolution, particularly that 
of a chaotic system. In a seminal paper, Peres 0] no- 
ticed that quantum time-evolution can be sensitive to 
the differences between chaotic and integrable dynamics 
in a peculiar set up: One needs to examine the overlap 
of identically prepared states, but evolved with slightly 
different Hamiltonians. This overlap, called Loschmidt 
echo (LE) or fidelity Q, measures the recovery obtained 
when a wave packet evolves for a time t, followed by a 
backwards evolution with a perturbed Hamiltonian for 
the same time interval. 

A considerable number of investigations has been de- 
voted to study the interesting and intricate phenomena 
related to the LE, in particular the different regimes that 
arise depending on the perturbation strength. For very 
small perturbations, the LE is described by standard 
perturbation theory and a Gaussian decay is observed. 
For stronger perturbations, where perturbation theory 
breaks down, large phase fluctuations Q lead to an ex- 
ponential decay of the LE described by the Fermi Golden 
Rule (FGR) 0O| For even stronger perturbations, but 
still weak in the classical sense, a semiclassical analysis 
yields an exponential LE decay that does not depend on 
the perturbation strength: The decay rate is determined 
by the Lyapunov exponent that characterizes the classi- 
cal counterpart of the unperturbed system Q . The latter 
two cases are called the FGR and Lyapunov regimes re- 
spectively, where the LE decay rate is the minimum be- 
tween the width of the local density of states (LDOS), as 
given by the FGR, and the Lyapunov exponent 0, 0, 13 • 
These predictions were verified numerically in a number 
of systems S S @- The theory is successful to 
the extend that, by analyzing the LE decay, the quan- 



tum evolution of a system can be used to quantitatively 
assess its classical Lyapunov exponent 0, |jj . 

The theory was later extended to classically integrable 
systems [Tol ] , in which case a power law like decay is pre- 
dicted. This result is still somewhat controversial [lT| . 
since as a rule integrable systems display non-generic fea- 
tures 0| . In any event, these works indicate that the 
LE decay is very different whether the underlying classi- 
cal system has a chaotic, integrable, or even mixed phase 
space [l|- 

Albeit this wealth of interesting results, so far the the- 
ory of the LE non-perturbative regime has only dealt 
with time independent perturbations. The most prob- 
able motivation for this restrictive choice can be traced 
back to the experiments that triggered the research on 
the LE problem 0, 0] : They studied the time reversal 
of many-spin dynamics, where the perturbation is simply 
a static part of the Hamiltonian. 

Numerous physical situations call for an extension of 
the LE theory that accounts for a time-dependent per- 
turbation. Let us explicit mention a few. Experimen- 
tally, a subsystem selected from a large spin system with 
many-body interactions can be represented as immersed 
in an external fluctuating potential [l6| - the same ap- 
proximation holds whenever the uncontrolled degrees of 
freedom are those of an environment with complex dy- 
namics. Formally, the current analytical description con- 
trasts with numerical results [ljobserved in periodically 
kicked one-dimensional models 0] , where the pertur- 
bation can be interpreted as time dependent. The need 
is further stressed by the relevance of the LE to quantum 
computation [Lsl llil |20| . decoherence in open systems 
and mesoscopic physics |24|. Indeed, the 
decay of the LE is related to the decay of quantum cor- 
relations and the quantum-classical correspondence, as 
can be shown using the Wigner function representation 

mm mug. 

In this work we use the semiclassical approximation to 
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derive the LE decay in the presence of a time-dependent 
perturbation, generalizing the approach presented in 
Ref . 13- We show that the existence of a LE perturbation- 
independent regime is quite generic. For that purpose, 
instead of using a particular model, we use a statisti- 
cal approach. We obtain a closed expression for the LE 
decay in the FGR regime using simple assumptions on 
the perturbation autocorrelation function. We conclude 
by discussing the different limits of our results and the 
seemingly strange feature that faster fluctuations of the 
perturbation or stronger chaos in the system lead to a 
slower decay of the Loschmidt echo 



II. LOSCHMIDT ECHO IN A 
TIME-DEPENDENT ENVIRONMENT 



The semiclassical propagator reads 



The object of interest, the Loschmidt echo, is defined 



as: 



M(t) = \(i; \U(to,t)U (t,toM W 



(1) 



where \ipo) is an arbitrary wave packet prepared at time 
to. For simplicity, and in line with Ref. 3, we choose 
the initial state \tpo) as a Gaussian wave-packet centered 
at an arbitrary point ro with dispersion a and initial 
momentum po. This restricted choice can be relaxed by 
considering other kinds of localized states in phase space 
|28L l2p . evolved states ^(Jj an d even eigenstates of Ho 
0, |26(. In Eq. U$ is the standard time evolution 
operator, namely 



U (t,t ) = Texp 



dt'H (t') , 



where T is the time ordering operator, while 

rio 



U(t ,t) = Texp 



dt'H(t'] 



(2) 



(3) 



with T the inverse time ordering operator. Equation 
is also viewed as the fidelity of two wave packets prepared 
at the same initial state and evolving forward in time 
under different Hamilton operators. 

In general, time ordering makes the exact evaluation of 
M(t) for a time dependent Hamiltonian a daunting task. 
To circumvent this difficulty we employ the semiclassi- 
cal approximation, in which time ordering is trivially 
accounted for by taking the time evolution of classical 
trajectories, as we detail in the sequel. 

We consider the Hamiltonian H defined as 



H = H + V(q,t), 



(4) 



where Hq is a time independent Hamiltonian that dis- 
plays chaotic motion in the classical limit and T^(q, t) is 
the time-dependent perturbation potential or the system 
interaction with a complex environment. 



(q'l^(*)|q) 



2nhi 



d/2 



s(q',q,t) 



xexp ( -5 s (q',q,t)- ^a s ) , (5) 



where s is a classical path that spends a time t to travel 
from q to q', S s is the action (Hamilton principal func- 
tion), given by fi^q', q,i) = J Q dri(q s (r), q s (r), r), a s 
is the number of conjugate points along s, and C s is 
the Jacobian of the phase-space transformation between 
<5p'(0) and <5q'(i) - a density of classical paths. 

It is only possible to proceed analytically if we re- 
strict ourselves to the regime of weak perturbations, in 
the sense that classical perturbation theory is applica- 
ble. More specifically, we approximate the action along 
a given trajectory s by 



S s (t) ~ S° g (t) + [ dt'V(<i s {t'),t'), 
Jo 



(0) 



where S°(t) refers to the action corresponding to s ob- 
tained from H and q s (t) gives the particle position along 
the unperturbed trajectory s as a function of time. For 
chaotic systems, this approximation is accurate up to a 
time t cp proportional to the logarithm of the strength 
of V . In this sense, the perturbation is weak when t cp 
becomes the largest time scale of the problem. This re- 
striction does not preclude the perturbation to be quan- 
tum mechanically large, since the actions are measured in 
units of h [5j • It has been observed that the classical per- 
turbation approximation works, in general, surprisingly 
well even for times longer than t cp . This has been related 
to the structural stability of the manifold of trajectories 
in phase space 30]: Even though individual trajectories 
are exponentially sensitive to perturbations, one can al- 
ways find a "replacement" trajectory in the manifold that 
joins the points of interest for a given time interval [231. 

Our calculation proceeds along the lines of Ref. y, 
which we now briefly sketch. We assume that the wave 
packet (r\ipo) is well localized, £ ^> a ^> XdB, where £ is 
a typical length of the perturbation (in Ref. y the width 
of Gaussian impurities) and XdB is the de Broglie wave- 
length of the particle. Neglecting terms with a rapidly 
oscillating phase, one arrives at the semiclassical expres- 
sion for the Loschmidt echo, 



x exp 



dr C s exp 

s(r,r ,t) 



n 



h 2 



(7) 



where AS S is the action difference between trajectories 
evolved with H and H, and p = — dS s /dr\ r = ra . All tra- 
jectories s start at Tq, the position where the Gaussian 
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wave packet (r|-0o) is centered at. For short times M{t) 
can only capture the local instabilities of the classical dy- 
namics, thus, it shows large fluctuations p|. By sampling 
over different initial values of i"o and po, or over an en- 
semble of perturbations, one obtains an average (M(t)) 
that puts in evidence the exponential decay |25| . The 
Lyapunov and the FGR decay regimes are related to the 
different ways of pairing the path summations in the dou- 
ble sum of Eq. J7J . 



A. Non diagonal contributions to (M(t)) 

Let us first calculate the terms where the two trajec- 
tories lie far apart in phase space. Such contributions to 
(M(t)) are usually called non-diagonal (different trajec- 
tories), and read 



(M nd (t)) 



^2 



x exp 



dr ^ C s 

s(r,r ,t) 



° I— \ 2 

Jfl (P.-Po) 



exp 



-AS s (t) 



(8) 



where (...) indicates that we average over the wave packet 
initial positions ro, as well as over an ensemble of pertur- 
bations. 

We assume, as is customary for chaotic systems, that 
the actions for different paths are uncorrelated and Gaus- 
sian distributed [28l l3l| . This leads to an enormous sim- 
plification, allowing us to write 



exp 



7M s (i) 

n 



exp 



1 



(9) 



We remain with the task of evaluating the action variance 



({As s (t)} 2 )= / at' at"(v(q s (t'),t')v(<i s (t")X)). 



(10) 

For that purpose we introduce an ensemble of perturba- 
tions V to model the general features of the environment. 
We replace the phase space average (• • • ) by the ensem- 
ble average 7TT , the equivalence between averages being 
supported by the ergodicity of the system. In order to 
keep our calculation as general as possible, we assume 
very little knowledge of the perturbation, requiring only 
that time and space correlations arc independent, viz. 



V(q, t)V(q', t>) = V 2 C s (|q - q!\)C T (\t - t'\ 



(11) 



The typical perturbation strength is (V 2 ) 1 ' 2 , and ry = 
hjiy 2 ) 1 ! 2 is its associated time scale. The dimensionless 
functions Cs and Ct quantify the spatial and time corre- 
lations of the potential V(q, t). We further require that 
Cs or Ct decay sufficiently fast, so that 



Arr d - l C s {r) < oo 



and / dtC T (t) < oo. (12) 
Jo 



For chaotic systems this is a sensible assumption. 

To guide the discussion, let us introduce the correla- 
tion length £o and the correlation time To that charac- 
terize Cs and Ct respectively. Since the average (tTTJl) is 
computed along the classical trajectories of the system, 
the asymptotic decay l|12l) can be induced not only by 
the fluctuations of V, but also by the intrinsic chaotic 
dynamics of Hq. In general, £o and ro are given by the 
minimum between the natural scales of V and Ho. For 
instance, when the perturbation is a static change in the 
mass tensor of a free particle bouncing off the walls of 
a billiard system, £o is solely given by the dynamics of 
Hq and is equal to the mean free path between collisions 
IE 0H3 Another example can be found in Refs. 
and|& where the effective scale ro is given by the kicking 
period of the unperturbed Hamiltonian - although the 
perturbation is a time independent change in the kicking 
strength. Hence, our results are valid not only for ran- 
dom perturbations, but also for static and periodic ones: 
the chaoticity of the underlying Hamiltonian alone can 
enforce conditions l|12|) . 

In the limit of ro S> 1/A the perturbation is quasi-static 
and the results of Ref. 3 hold without further change. We 
are interested in the regime where the typical times of 
the perturbation are comparable to those of the system, 
ro < 1/A. 

Replacing space averages by ensemble averages (|TTj> . 
we write Eq. I jlOj l as 



(AS s (t) 2 ) 



V 2 



x C R 



dt 



dr 



OO 
T „ 



<ls(i-~)-q s (t+-) 



C t (t$3) 



where we considered times t much larger than ro and 
£o/v, which allows us to take the integral in r from — oo 
to +oo. Eq. (|13fl has two limiting regimes that are readily 
solved. In the first one, the spatial disorder has a much 
shorter scale than the temporal one: ro ^> £o/ u = r^ . In 
this case the decay of M nd (t) is dominated by the same 
exponent as the one found in Ref. |jj 



(AS s (t) 2 ) 



V 2 f dt 
Jo 



drC s 



q s (*-^)-q s (* + ^) 



(14) 



where Ct(t) is assumed constant and f is given by a 
FGR calculation 

1 

fi 



(15) 



When ro <C r^, we deal the opposite regime, and 



(AS s (t) 2 ) ~V 2 dt 
Jo 



drC* T (r) = — K 2 

T2 



with 



1 

r 2 



T 2 ' 



(16) 



(17) 
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Thus, in these two limits and complementary situa- 
tions, the FGR exponent changes from being governed 
by the spatial to the temporal correlations of V(q,t). 
The interesting "correlation crossover regime" - where 
neither the temporal nor the spatial correlation domi- 
nate - will be discussed shortly for a particular form of 
Cs and Ct- 



B. Diagonal contributions to (M(t)) 

Let us first explicitly write (JJJ, namely 



M(t) 



(r.iD.t) 
'(r', r ,t) 



x yexp 
x exp 



(AS s (t) - AS s ,(t)) 



2 

-^2 ((Ps - Po) 2 + (P s - Po) 



and analyze the case where the trajectories s and s' re- 
main close to each other. Now the action differences can- 
not be considered as uncorrelated, and we have to take 
into account the fluctuations in 



exp 



(AS a (t)- AS A*)) 



exp 



2h 2 



-AS s (t)-ASs>(t)Y 



In the same order of approximation of Eq. (jSJ, we write 

AS s (t)-AS s >(t)= [ df [V(q s (t'),t')-V(q s ,(t'),t')}. 
Jo 

(20) 

As the two trajectories remain close in coordinate space, 
we can expand V(q s (t),t) to first order around s and 
obtain 

AS s (t)-AS s ,(t)~ f dt'VV(q s (t'),t')-[q s (t')-q s ,(t')] 
Jo 

To calculate the action difference variance we turn our 
attention to the force correlation function, namely 

Cv(|q ~ q'|, \t - t'\) ee (W [q, t] ■ W [q', t'}) . (22) 

As before, we introduce an ensemble of perturbations, 
and write 

Cv(|q-q'|, \t-t'\) = V ff C r (|t-*'|)(V q -V q /)C lS (|q-q , |), 

(23) 

such that (V q • V q <)Cs (|q — q'|) decays sufficiently fast, 
in the sense defined by Eq. \V1\- 

As time evolves, the separation between the coordi- 
nates q s (i) and q s '(t) grows as e At , where A is the largest 



Lyapunov exponent of Ho. As a result, after some alge- 
bra, Eq. (|T^|l gives exp [— A(r — r') 2 /h 2 }, with 



A = V 2 



TO- 



-2\t 



2A 



(24) 



when Ct dominates the decay of Cy , and 



A = V 2 



1 



-2At 



2\v 



dq 



l-ddC's(q) d 2 C S (q) 



dq 



dq 2 



(25) 

when Ct decays slowly. 

In summary, the main result of Ref. holds, namely 

M(t) =Aexp(-Xt) + B exp(-t/T), (26) 

where A = [ma/(A 1 / 2 t)] d , A is the classical Lyapunov ex- 
ponent of the system and 1/f is given by Eq. I|13|l . The 
exponential decay of the LE is dominated by the small- 
est between 1/f and A, giving a crossover from FGR to 
Lyapunov decay as the perturbation strength increases. 



C. Correlation crossover 

In the regime where r w t£, one can only obtain fur- 
ther insight by assuming a specific form of the correlation 
functions. Although it is a less general result, one can still 
encompass a broad class of possible perturbations whose 
correlator decay in a particular way. We will consider the 
case where both Cs and Ct have Gaussian shapes, 



(19) V(ci,t)V(q',t') = 



V 2 



■ exp 



,'|2 



exp 



\t-t'\ 



t\2 



Under the assumption that t is large compared to To and 
T£, we replace in Eq. (|13fl and Eq. (|19fl . and obtain the 
decay rate for the FGR regime 



1 

f 



and the prefactor A of the Lyapunov regime: 



A = 



v 2 \r- 



1 -e 



-2At 



d-1 



T 2 T 2 



(28) 



(29) 



When the temporal or spatial correlation dominate, we 
recover the previous limit 



1 T c 

— ~ — with r c = minlr. 



0; '£J 



(30) 



As before, if the effective time scale r c becomes too short, 
the perturbation cancels itself out causing a very slow de- 
cay. This result is consistent with studies of time depen- 
dent errors in a quantum computer 20], where the dy- 
namical decoupling to the environment was interpreted 
as a manifestation of the quantum Zeno effect 0, . 
Notice that when r c is dominated by the dynamics of Hq, 
the fluctuations become faster for chaotic systems with a 
larger A (25|. 
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III. CONCLUSIONS 

We have extended the semiclassical theory of the 
Loschmidt echo to cope with time dependent perturba- 
tions. We expect our results to remain valid in more 
complex or analytically difficult cases, suitable only for 
numerical studies. Our treatment is sufficiently general 
as to describe the situations where the perturbation is 
the random effect of an uncontrolled environment on the 
system. The fluctuations we considered could arise either 
from an explicit time dependence of the perturbation po- 
tential, or from the ergodic nature of Hq- In the last case, 
the underlying chaotic dynamics mimics the randomness 
required for the decay of the correlation functions. Thus, 
our results should also apply to periodic or very simple 
oscillating perturbations. 

We showed that the Loschmidt echo Lyapunov regime 
is barely affected by the time-dependence of the pertur- 
bation, except for prefactors: The decay is dominated by 
the system's intrinsic dynamics of stretching and folding. 
In the FGR regime - when the non-diagonal terms dom- 
inate - the spatial and time scales of the perturbation 
compete with each other, and a simple behavior can be 
extracted when the relevant scales are far apart. In the 



intermediate regime, where the scales are comparable, us- 
ing a simple (yet general) example we compute the decay 
rate of M(t). The form of Eq. H3()|l stresses how fast fluc- 
tuations lead to self-cancellation of the interaction with 
the environment. In the case of the LE, a vanishing FGR 
exponent prevents the appearance of the perturbation 
independent Lyapunov regime. Surprisingly, this hap- 
pens not only for rapidly fluctuating perturbations, but 
also by increasing the Lyapunov exponent. The slowing 
down of the FGR regime of decoherence - induced by fast 
fluctuations - was recently experimentally measured in 
NMR experiments [3j| , where a connection to the quan- 
tum Zeno effect was observed. It is interesting to recall 
that dynamical decoupling to the environment is what 
makes liquid NMR quantum computers possible (albeit 
small). The fast random movements of the molecules in 
the liquid average out the more difficult to control dipo- 
lar interactions present, e.g., in solids. Our work points 
to the importance of exploring dynamical alternatives to 
suppress quantum decoherence |32(. 
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port from CNPq (Brazil) and CONICET, ANPCyT and 
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